Low temperature specific heat of ID multicomponent systems at commensurate- 
incommensurate phase transition point 
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Low temperature dependence of specific heat of one- dimensional multicomponent systems at the 
commensurate- incommensurate phase transition point is studied. It is found that for canonical 
systems, with a fixed total number of particles, low temperature specific heat linearly depends on 
temperature with a diverging prefactor. 
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Second order phase transitions in one- dimensional 
quantum or two- dimensional statistical systems arc often 
classified by fixed point Hamiltonians enjoying conformal 
invariancei. In such situations, low temperature thermo- 
dynamic behavior is fixed solely by the central charge^, c, 
the quantity that roughly measures numbers of degrees of 
freedom, and velocity of the linear excitation, v - 'speed of 
light'. Low temperature specific heat combines these two 
quantities and has the universal form^ C(T) = ircT/3v, 
depending linearly on temperature. The notable excep- 
tion from the above picture is provided by the second 
order phase transition driven by chemical potential cou- 
pled to the conserved quantity^, so called commensurate- 
incommensurate (C-IC) phase transition^. At the C-IC 
phase transition point system, in contrast to conformal 
invariance, obeys zero scale- factor universality-^ with dy- 
namical critical exponent z = 2. In particular for the spe- 
cific heat it implies square root temperature dependence 
for T -> 0. 

Describing the C-IC phase transitions in multicom- 
ponent systems (systems of fcrmions, bosons or spin 
S chains) one can not rely on mode- decoupling 
approximation^. In particular, as a result of spin- charge 
non- separation, magnetic susceptibility of spin gapped 
systems, at the edge point of magnetization plateau, 
stays finite instead of diverging unless special microscopic 
symmetries are present (e.g. particle-hole symmetry at 
half filling for attracting Fermi Hubbard model). This 
kind of behavior has also been observed for a number of 
integrablc models such as the two component Fermi Hub- 
bard model both on lattic e) 9 ' 10 as well as in continuum^, 
or spin S generalization of the integrable t — J chain 
doped with 5—1/2 carriers^. Similar behavior holds for 
generic systems of multicomponent Fermi or Bose gases, 
with gapped spin and gapless charge excitation&i^. No- 
table example of bosonic systems, in this respect, is pro- 
vided by antiferromagnetically interacting spin 1 Bose 
gas^ii. 

Zero temperature magnetization increases linearly 
with the field at the edge of the magnetization plateau, 
due to the finite value of magnetic susceptibility. This 
kind of a behavior is characteristic to conformally invari- 
ant systems. The natural question arises whether some 
kind of conformally invariant theory can be applied to 
describe the edge points of magnetization plateaus. For 



this reason we study temperature dependence of specific 
heat in multicomponent systems at C-IC phase transition 
point where modes do not decouple. We find that if total 
number of particles is fixed then square root temperature 
dependence, which holds if one works with fixed chem- 
ical potential and also is characteristic to mode- decou- 
pled systems undergoing C-IC transition (for both cases 
with fixed chemical potential or fixed total number of 
particles, e.g. attractive Hubbard model at half filling 
at the onset of magnetization) , is modified by linear law, 
however the prefactor diverges as In 2 T for T -> 0. To 
our knowledge no previous studies on finite temperature 
thermodynamic behavior at C-IC phase transition point 
have been reported for multicomponent systems where 
modes do not decouple. 

For concrete calculations we will consider a system of 
two component attractive fermions and place it in an 
external magnetic field equal to the value of the spin 
gap. However our main claim on the disappearance of 
the square root temperature dependence of the specific 
heat (at C-IC phase transition point) will hold true for 
generic situation, where gapped mode (spin) , undergoing 
non- rclativistic softening at the critical point, couples to 
linearly dispersing excitations (charge). This statement 
will apply to all situations where magnetic susceptibility 
stays finite at the edge points of magnetization plateaus 
in zero temperature limits. 

To simplify things drastically we will consider the di- 
lute limit of the system of attractively interacting two 
component fermions. Since the dilute limit falls under 
strong coupling, the ground state (for T = 0) is made 
of bound pairs and the low temperature thermodynamic 
properties can be modeled by the mixture of noninteract- 
ing pairs and thermally created un-compensated spin- up 
particles (we choose finite temperature magnetization to 
be positive for nonzero magnetic field). Spin- down par- 
ticles will also be created thermally, but since they have 
spectral gap (which is large in strong coupling) their den- 
sity will be exponentially suppressed at low temperatures 
and thus they will be ignored in the following. For in- 
tegrable Fermi Hubbard model (both on lattice^ and 
in continuum) the thermodynamic Bethc ansatz method 
can be applied^, and in particular recently it was shown 
that in dilute limit this method simplifies considerably^. 
However, for finding specific heat at the critical point 
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when magnetization sets in we need not use intcgrability 
so that our reasoning will stay both simple and general. 

Phcnomcnological Hamiltonian describing the low en- 
ergy properties of the two component attractive fermions 
in strong coupling at the C-IC critical point reads: 



H = J2(E P (k) - E F )a\a k + £ E^k)c{ 



Ck, 



(1) 



where a){c)) and a(c) are pair (up- spin fermion) cre- 
ation and annihilation operators. For calculating ther- 
modynamic properties it does not matter whether pairs 
are modeled by hard core bosons (which they are) or 
fermions^. 

At the edge of the zero magnetization plateau disper- 
sion of up- spin particles is approximated by a purely 
quadratic expression: 



E f (k) = yj- 



v*k 2 



A 



w|fc72A. 



(2) 



Here A is spin gap opened due to strong attraction be- 
tween spin- up and spin- down fermions, and u-f is the 
velocity of gapped excitations which, as we will see later, 
will not appear in the leading temperature dependence of 
specific heat. Since low temperature thermodynamics is 
determined by low lying modes (states near Fermi wave- 
vector) for the dispersion relation of pairs we can use a 
simplified strictly linear expression: 



E p (k)-E F = v p (\k\-k F ). 



(3) 



Fermi wave vector k F is related to the linear density of 
pairs n p by the usual ID relation k F — n p Tr. In the follow- 
ing n p will stand for the density of pairs strictly at zero 
temperature: n p = n p (T = 0). Dispersions E^(k) and 
Ep(k) — E F considered in our simplified model are plot- 
ted on Fig.([T]). At zero temperature the ground state is 
occupied solely by pairs, and density of un-compensated 
up- spin electrons is zero. Although one might argue 
that Hamiltonian ([T]) lacks microscopic derivation, it is 
expected to exactly reproduce leading low temperature 
dependence of the specific heat at C-IC phase transi- 
tion point in dilute limit, the same way as it yields exact 
expression for magnetic susceptibility^. The thermody- 
namic potential of the grand canonical system can be 
obtained straightforwardly: 
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For specific heat one obtains easily that the leading tem- 
perature dependence is linear from pairs (usual confor- 
mal contribution^ with central charge c = 1 and veloc- 
ity of pair motion v p ) and square root (which dominates 
over the linear one in the low temperature limit), from 
the up- spin electrons, due to the quadratic low energy 
dispersion they enjoy. This kind of behavior is also re- 
alized for canonical systems at special fillings, e.g. the 
attractive Fermi Hubbard model at half filling, where 
spin and charge modes decouple. At half filling charge 




FIG. 1: simplified dispersions with purely quadratic and lin- 
ear behavior adopted in this paper. If total number of par- 
ticles is conserved creation of each pair of up- spin fermions 
(smaller bullets) leaves a hole in the Fermi sea of bound pairs 
(larger bullet). 



sector of Hubbard model is described by conformally in- 
variant SU\{2) Wess- Zumino models, while the spin 
sector, when magnetic field equals to singlet- triplet spin 
gap, is characterized by quadratic low energy behavior 
responsible for the leading, square root, temperature de- 
pendence of specific heat. In the grand canonical systems 
at finite temperatures, the thermal fluctuations create a 
finite density of up-spin electrons given by the following 
expression: 

n t (T) = -C(l/2)V2vrA(v / 2 - l)VT/2nv^, (4) 

where C(l/2) is the value of the Ricmann £ function at 
0.5. The number of pairs also gets modified by tempera- 
ture fluctuations accordingly: 



h p (T) 



T 



TTV V 



k F /T< 



.(5) 



Thus, due to finite ground state population of pairs, ther- 
mal fluctuation of density of pairs is exponentially sup- 
pressed, whereas for up- spin electrons, due to Van Hove 
singularity in their density of states, square root depen- 
dence holds. Already from the above picture it is clear 
that working at a fixed total number of particles can mod- 
ify the leading square root behavior of specific heat at C- 
IC point, since thermal occupation of up- spin band will 
be constrained by the reluctance of depletion of the band 
of bound pairs. Since each broken pair produces two up- 
spin electrons, as depicted on Fig.([T|), the constraint that 
fixes the total number of particles reads: 



n t (T) + 2n p (T) = 2n p . 



(6) 



To implement the constraint we introduce the temper- 
ature dependent chemical potential, which we determine 
from the equation: 



dk/2-K 



1 



dk/n 



1 



2n p . (7) 



The consistent solution of the above equation requires 
that jj,(T)/T — > — oo, when T — > 0, so that we obtain: 



M (T) ~ -aVTe^/ T , 



(8) 
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with positive coefficient a 



In the zero tem- 



perature limit we obtain the following solution for the 
chemical potential: 
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Now we will use the obtained chemical potential Eq.© 
to calculate the grand canonical potential: 
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Dots stand for sub-leading terms at low temperatures. 
Grand canonical potential is related to energy, E, and 
entropy S, by the following thermodynamic relation: 

il = E -TS - n(T)N t (T) - (v p k F + 2n(T))N p (T). (11) 

For the leading temperature dependence of energy we get 
the following expression: 



E 



(v p k F + 2/i(T)) 2 irT 2 27>(T) 
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For thermally depleted density of pairs we obtain: 

frr s v P k F + 2/i(T) 
n p\ T ) = — > 



(12) 



(13) 



and the density of thermally created up- spin electrons is 
given by the following expression: 



n t (T) = 2n p - 2n p (T). 



(14) 



Note, due to constraint of keeping total number of parti- 
cles fixed, the square root dependence Eq.Q is modified 
in the following way: 



n t (T) 



2Tln(T/A) 



T/A < 1. 



(15) 



Using thermodynamic relation Eq. ([lip we obtain for 
the leading temperature dependence of entropy the fol- 
lowing expression: 



s(T) = S(T)/L 



4/i 2 (T) 6/i(T) ttT 



Tttv v 
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(16) 



Finally substituting chemical potential from Eq. © for 
the leading temperature dependence of specific heat we 
get: 



C(T) = T 



ds(T) _ Tin 2 (T/A) 
dT irv n 



(17) 



in the limit T — > 0. This is our main result: specific heat 
depends linearly on temperature, albeit with a diverg- 
ing prefactor. This is kind of an intermediate between 
the square root and linear dependences. Recent stud- 
ies reported similar linear behavior of the specific heat 
with diverging prefactor for critical two- subband quan- 
tum wirei^ when the second band started to fill. How- 
ever in our case we can not interpret this as a signal of 
restoration of Lorentz or conformal invariance. 

At this point we would like to comment on the differ- 
ences between working at fixed chemical potential and 
fixed total number of particles for attractive Hubbard 
model away of half filling at the onset of magnetization 
(C-IC transition point). For the case of fixed total num- 
ber of particles magnetic susceptibility stays finite and, as 
we showed here, specific heat depends linearly on tem- 
perature (with diverging prefactor). However, for the 
case of fixed chemical potential magnetic susceptibility 
diverges and specific heat shows square root tempera- 
ture dependence- characteristic to the mode decoupled 
theory. One may suspect that if one works with the fixed 
chemical potential one can use mode decoupled theory, 
result of linear bosonization, to describe the C-IC phase 
transition point. For calculating static magnetic suscep- 
tibility and specific heat this statement is true. However, 
using mode- decoupled theory across the C-IC transition 
is misleading, because away of half filling and for non-zero 
magnetization spin and charge modes couple independent 
of whether one fixes chemical potential or total number 
of particles. In particular spin and charge stay coupled 
(spin- charge mixing angle stays finite) even in the limit 
of vanishing magnetization (C-IC point) for both cases 
with fixed total number of particles or fixed chemical po- 
tential. Within spin- charge decoupled theory one can 
not calculate correctly e.g. assymptotics of static sin- 
gle particle correlation functions^ in the limit of vanish- 
ing magnetization (C-IC point). Fixing total number of 
particles further enhances mode coupling and leaves its 
traces on magnetic susceptibility as well as specific heat 
(topic of this manuscript). 

The presented analysis is directly applicable also to 
the phase transition between partially and fully polarized 
states, when the down- spin fermions band gets empty 
(at h = h sa t, where h sat is magnetic field needed to com- 
pletely polarize the system). In this case one simply has 
to replace constraint Eq.© by: N^(T) + N±(T) = const, 
thus specific heat at h = h sat will have the similar tem- 
perature dependence as in Eq . ([!?)) . 

For completeness we write out the leading low tem- 
perature dependence of specific heat in two component 
attractively interacting fermions in an external magnetic 
field for all values of the field away from half filling: 

{T/v c h < 2A or h > h sat 

Tin 2 (T/A) h = 2A or h = h sat (18) 
T(v c + v s )/v c v s 2A < h < h S at 

where v c and v s denote charge and spin velocities^. Note 
that v s (m) ~ m when m — > 0, and v s (m) ~ n — 2m when 
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m — > n/2, where m is the magnetization density and n 
the total density of particles. In the region A < h < 
h sa t spin and charge modes are coupled, and they do 
not get separated even in the limits h — > 2A + and h — > 
h sat _£*m. This non- separation of spin and charge modes 
explains disappearance of the square root temperature 
dependence of specific heat at critical points h = 2A and 
h = h sai r^. Exactly at half filling, where spin- charge 
separation holds, C(T) ~ s/T for h = 2 A or h — h sat . 

To summarize on simple example of two component 
attractive fermions placed in an external magnetic field 
equal to the spin gap we studied the thermodynamic 
properties of C-IC phase transition point in multicom- 
ponent systems. Our findings are generic, namely, we 
expect that for multicomponent systems square root tem- 
perature dependence, which holds for the fixed chem- 
ical potential, will be changed by linear dependence 
(with diverging prefactor) if instead total number of 
particles is fixed. In particular, this statement equally 
applies to the edge points of magnetization plateaus 
of ID bosonic systems, i.e. spin 1 bosons interacting 
antifcrromagnetically^i. It also holds for the edge points 
of mid magnetization plateaus of multicomponent sys- 
tems, i.e. in doped dimcrized Hubbard chains^ 2 -. How- 
ever, currently we can not determine if the In 2 T prefactor 
in the expression of specific heat Eq. (fTT|) is universal. To 
answer this question, it will be useful to study e.g. the 



same system of two component attractive fermions in the 
weak- coupling limit. One can make use of integrabil- 
ity of Fermi Hubbard model there. Especially intriguing 
will be to look close to half filling to follow how does 
the square root dependence of specific heat cross over to 
linear behavior with diverging prefactor as a function of 
doping. 

Our findings may be relevant for real quasi one- 
dimensional spin gapped materials, where one can study 
experimentally specific heat as a function of applied mag- 
netic field and doping. Square root dependence, which 
is expected at half filling, should be modified by linear 
with logarithmic prefactor, away from half filling, at crit- 
ical magnetic field strength equal to the spin gap. One 
can as well measure the leading temperature dependence 
of specific heat at saturation magnetization in a generic 
ID electron system away of half filling and compare the 
results to Eq. ((T7|) . 
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